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ABSTRACT 

We report multi-epoch VLBI observations of the source 1741—038 (OT— 068) as it 
underwent an extreme scattering event. Observations at four epochs were obtained, 
and images were produced at three of these. One of these three epochs was when the 
source was near the minimum flux density of the ESE, the other two were as the flux 
density of the source was returning to its nominal value. The fourth epoch was at the 
maximum flux density during the egress from the ESE, but the VLBI observations had 
too few stations to produce an image. 

During the event the source consisted of a dominant, compact component, 
essentially identical to the structure seen outside the event. However, the source's 
diameter increased slightly at 13 cm, from near 0.6 mas outside of the ESE to near 
1 mas during the ESE. An increase in the source's diameter is inconsistent with a 
simple refractive model in which a smooth refractive lens drifted across the line of 
sight to 1741—038. We also see no evidence for ESE-induced substructure within the 
source or the formation of multiple images, as would occur in a strongly refractive 
lens. However, a model in which the decrease in flux density during the ESE occurs 
solely because of stochastic broadening within the lens requires a larger broadening 
diameter during the event than is observed. Thus, the ESE toward 1741—038 involved 
both stochastic broadening and refractive defocussing within the lens. If the structure 
responsible for the ESE has a size of order 1 AU, the level of scattering within an ESE 



lens may be a factor of lO'^ larger than that in the ambient medium. A filamentary 
structure could reduce the difference between the strength of scattering in the lens and 
ambient medium, but there is no evidence for a refractively-induced elongation of the 
source. We conclude that, if ESEs arise from filamentary structures, they occur when 
the filamentary structures are seen lengthwise. 

We are able to predict the amount of pulse broadening that would result from a 
comparable lens passing in front of a pulsar. The pulse broadening would be no more 
than 1.1 fis, consistent with the lack of pulse broadening detected during ESEs toward 
the pulsars PSR B1937+21 and PSR J1643-1224. 

The line of sight toward 1741—038 is consistent with a turbulent origin for the 
structures responsible for ESEs. The source 1741—038 lies near the radio Loop I and 
is seen through a local minimum in 100^ emission. 

Subject headings: ISM: general — quasars: individual (1741—038) — radio continuum: 
ISM — scattering 



1. Introduction 

Extreme scattering events (ESE) are a class of dramatic changes in the flux density of radio 
sources (Fielder et al. 1994a). They are typically marked by a decrease (^ 50%) in the flux 
density near 1 GHz for a period of several weeks to months, bracketed by substantial increases, 
viz. Fig. 1. Because of the simultaneity of the events at different wavelengths, the time scales 
of the events, and light travel time arguments, ESEs are likely due to strong scattering by the 
Galactic interstellar medium (ISM; Fiedler et al. 1987a; Romani, Blandford, & Cordes 1987). 
First identified in the light curves of extragalactic sources, ESEs have since been observed during 
a timing program of the pulsars PSR B1937+21 (Cognard et al. 1993; Lestrade, Rickett, &: 
Cognard 1998) and PSR J1643-1224 (Maitia, Lestrade, k Cognard 1998). 

Extreme scattering events are generally ascribed to be the result of strong interstellar 
refraction by discrete ionized structures (Fiedler et al. 1987a; Romani et al. 1987; Clegg, 
Fey, & Lazio 1998), though Fiedler et al. (1994a) developed a model in which ESEs are due 
to stochastic broadening of source diameters. The extent to which refractive defocussing or 
stochastic broadening dominates during an ESE may be determined through VLBI imaging of a 
source undergoing an ESE. If refractive defocussing dominates, the position of the source may 
wander, and the shape of the lens or substructure within the lens may produce distortions in the 
appearance of the source; if the refraction is strong enough, the source may be multiply imaged 
(Cordes, Pidwerbetsky, & Lovelace 1987; Clegg et al. 1998). If stochastic broadening dominates, 
the source structure should remain largely unchanged except for an increase in its diameter. The 
extent to which one or the other of these mechanisms dominates may also provide clues to the 
origin of these lenses. To date, the only observational constraints on the ESE mechanism — besides 



the light curves — have been the lack of pulse broadening (as might be expected from a stochastic 
broadening model) and a variation in the pulse time of arrival (as would be expected from a 
refractive defocussing model) during the ESEs toward PSR B1937+21 (Cognard et al. 1993; 
Lestrade et al. 1998) and PSR J1643-1224 (Maitia et al. 1998). Of course, it need not be the case 
that the ESE-like events observed toward these pulsars resulted from the same kind of structure 
responsible for the ESEs toward extragalactic sources. 

This paper reports the first VLBI observations of a source (1741—038, OT— 068) while it was 
undergoing an ESE. In §2 we describe the observations, in §3 we discuss the implications of our 
observations, and in §4 we present our conclusions and suggestions for future work. 



2. Observations 

Figure 1 shows a portion of the 2.25 GHz light curve of 1741—038 as obtained by the US 
Navy's extragalactic source monitoring program at the Green Bank Interferometer (Fiedler et 
al. 1987b; Waltman et al. 1991; Waltman et al. 1999). Clearly evident is an approximately 
50% decrease in the source's flux density. The minimum occurred on or near 1992 May 25 
(JD 2448768.264), and the ESE is nearly symmetric about this epoch. We show only a portion 
of the light curve in order to focus on the ESE. The complete GBI light curve of 1741—038, 
extending from 1983 to 1994, has been published previously (Clegg et al. 1998). 

The epochs of observations are indicated in Figure 1, and the observing log is given in Table 1. 
Also shown are the GBI flux density measurements closest in time to the VLBI observations. 
While the GBI measurements were not simultaneous with the VLBI observations, the largest 
elapsed time between the GBI and VLBI observations was no more than 1 day. All observations 
were conducted in VLBI Mk II mode and recorded left circular polarization only. As part of this 
program, observations were also obtained on 1993 June 2, well after the ESE concluded. 

Shen et al. (1997) conducted 6 cm VLBI observations in 1992 November. Fey, Clegg, &: 
Fiedler (1996a; see also Fey, Clegg, &: Fomalont 1996b; Fey & Chariot 1997) have already presented 
images of 1741—038 at 6 cm obtained on 1989 April 12 and at 3.6 and 13 cm obtained on 1994 
July 8. In all images, the source is dominated by a single compact component, and there is little 
evidence of structural change within the source during the interval 1991-1994. These images have 
a higher dynamic range than we could obtain from our 1993 June observations. Consequently, we 
do not include the 1993 June observations here, but instead use the images from the other groups 
to discuss source characteristics outside the ESE. 

Nominal system temperatures for 1992 were obtained from the various stations. These were 
used for the initial amplitude calibration. A refined amplitude calibration was then determined in 
the following manner. The source 2121+053 (OX036) is a compact source used as a fringe finder 
during the observing programs. We fit a single, circular Gaussian to observations of 2121+053 



obtained in 1994 July/ finding a diameter of 0.9 mas. With this source structure and the source 
flux density measured from the GBI monitoring program, we then fit the 2121+053 data at each 
epoch for a single, antenna-based scaling factor. This antenna-based scaling factor accounts for 
the difference between the assumed and actual system temperatures. Typical corrections were 
10-20%. 

A crucial assumption, and a potential systematic error, of this method is that the source 
diameter of 2121+053 remained essentially constant over the interval 1992-1994. VLBI images 
of 2121+053 obtained at epochs bracketing the time of our experiment show it to be extremely 
compact (Wehrle et al. 1992; RRFID^). The GBI monitoring program also shows no significant 
flux density changes, such as a large increase in the flux density, indicative of structural changes, 
such as the emergence of a new component. As measured by the GBI, the flux density of 2121+053 
decreased by about 25% during 1992, from about 2 Jy to 1.5 Jy, as part of a longer term decrease 
following an outburst in 1988-1989. During 1993 the flux density stopped decreasing, remaining 
at about 1.5 Jy. Superposed upon those longer trends are shorter time scale variability; in the 
interval 1992 June-August, the flux density varied from near 2 Jy to as low as 1.3 Jy, then 
recovered to 1.7 Jy. 

Fringe fitting at 13 cm was performed in a two steps. The arrays used consisted of a 
reasonably close cluster of antennas (mostly in the southwest US) combined with a smaller number 
of far-flung antennas. We fit first for the fringe delays and rates of the antennas comprising the 
cluster, without solving for the delays and rates of the distant antennas. After applying these 
delay and rate solutions, we fit for the delays and rates for the distant antennas only but used all 
antennas in the fit. This fringe fitting procedure increases the probability that we will find fringes 
to all stations by focussing first on the stations with the highest signal to noise and those that 
should have the smallest rates and delays. 

After fringe fitting, a series of first, phase-only and then, amplitude-and-phase self-calibration 
iterations were used to account for short time-scale fluctuations in antenna gain amplitudes and 
phases. Because of the extremely compact structure of 1741—038 seen outside the ESE, we often 
used a point source as the model during phase self-calibration. Since the length of the GBI baseline 
(2.4 km) is much shorter than the shortest baseline in our VLBI array (the VLBA_PT-VLBA_LA 
baseline at 237 km), we used the GBI flux densities as zero-spacing flux densities when imaging 
the source. 

Unfortunately, because of either station or correlator problems during the 1992 August 6 
observations, only one hour of the 10-hour run had more than three stations on source 
simultaneously. With only three stations available for most of the run, we can do little more 
than phase-only self-calibration and fit simple models to the data. We shall therefore restrict the 



^ Observations obtained from the Radio Reference Frame Image Database (RRFID), 
(URL:http://www. usno.navy.mil/RRFID). 



majority of our comments to the images from the epochs 1992 June 8, June 20, and July 9. 

The amount of time on-source ranged between 6 and 10 . The resulting thermal noise limit 
on our maps is therefore approximately 0.5 mJy beam^^. The actual off-source rms noise levels 
are in the range 1-3 mJy beam^^ and are listed in Table 1. 

Figures 2~4 show the images for the epochs 1992 June 8, June 20, and July 9. We have fit 
one or more Gaussian components to the u-v data; the solutions to these fits are in Table 2. We 
discuss each epoch briefly and separately. 



2.1. 1992 June 8 (Fig. 2) 

At this epoch the source is near the minimum flux density. The source structure at this epoch 
is consistent with that seen outside of the ESE (1992 November, Shen et al. 1997; 1989 April, 
1994 July, Fey et al. 1996a), namely a dominant compact component and a weaker component 
to the south. We had difficulty finding a model that fit the data with a non-zero major axis for 
the secondary component, whereas Fey et al. (1996a) found a major axis of approximately 6 mas. 
This discrepancy could result from two effects. First, there could have been modest evolution of 
the secondary component between 1992 and 1994. Such evolution might also account for a portion 
of the modest source variability seen in the GBI monitoring of 1741—038. Second, our u-v plane 
coverage is not as extensive as that of Fey et al. (1996a). Consequently, our coverage may be 
sufficient to indicate the component's presence without allowing us to fit a detailed model to it. 



2.2. 1992 June 20 (Fig. 3) 

The source structure at this epoch consists of a single component. The secondary component 
is not apparent. Its apparent disappearance is the combination of two effects. First, the 
observation at this epoch has only four stations, so that detailed structure is likely to be lost. 
Second, the secondary component has an inverted spectrum, a ^ 2 {S^ oc z/". Fey et al. 1996a). 
Extrapolating this spectrum to 18 cm, we expect that the flux density of the secondary component 
will be 5-10 mJy, sufficiently weak that we would not detect it. 

In contrast to the 13 cm structure, the axial ratio at 18 cm is h/a = 0.1. We attach little 
significance to this low value, however. Because the array contains Hartebeesthoek, it is extended 
significantly in the north-south direction. 



2.3. 1992 July 9 (Fig. 4) 

At this epocii tlie flux density of the source is starting to return to its nominal value, but this 
epoch is prior to the peak flux density during the egress of the ESE. During the ingress and egress 
from an ESE, when the flux density is above its nominal value, ESE-induced changes in the source 
structure are most likely to be visible (Clegg et al. 1998). Since we are unable to image the source 
at its peak flux density (§2.4), this epoch represents our best chance for seeing any ESE-induced 
changes in the source structure. 

The source structure is little changed from that outside the event. The source continues to 
be dominated by a single compact component. The secondary component is not apparent in the 
image at this epoch. The image at this epoch has the highest ofF-source noise level of the three 
epochs. This high noise level is the result of numerous intervals during the observing run in which 
fringes could not be found. Thus, this observation was more like a series of "snapshots" rather 
than a pointed observation. The combination of high noise level and poor u-v coverage could 
contribute to difficulty in detecting the secondary component. 

Although limited u-v coverage and the high noise level in this image are at least partially 
responsible for the absence of the secondary component, another possibility is that the ESE lens 
did not cover the source fully. Just such a possibility is indicated from Clegg et al.'s (1998) results. 
They found that a lens comparable in diameter to the compact component, 0.5 mas, produced the 
best match to the light curve. The relatively simple structure of 1741—038 does not allow us to 
place any constraints on the axial ratio of the ESE lens, though. 



2.4. 1992 August 6 

As noted above, for the majority of this epoch, only three stations were observing the 
source. Model fitting to the available data indicates that the source continues to have a compact 
component. The model fitted to the visibility data is listed in Table 2. 



3. Discussion 

In this section we use our VLBI images of 1741—038 to infer various properties of the lens 
responsible for this ESE, and in particular, the extent to which refractive and difFractive scattering 
were important. We begin by showing that the difFractive scattering, as manifested by additional 
angular broadening oF 1741—038 during the event, did occur. Clegg et al. (1998) have modeled 
the 1741—038 as due solely to retractive defocussing by an ionized cloud. We discuss the (limited) 
extent to which our observations can test the predictions of their model. We then turn to the 
question of the mechanism by which ESEs can occur — refractive defocussing (e.g., Romani et 
al. 1987; Clegg et al. 1998) or stochastic broadening (Fielder et al. 1994a). We conclude the section 



by considering what the shape of our images impUes about the shape of the lens and what the hne 
of sight to 1741—038 imphes about the genesis of the lens. 



3.1. DifTractive Properties of the 1741-038 ESE Lens 

Our observable for studying the diffractive effects of this lens is the angular broadening of 
the compact component of 1741—038. Plasma density fluctuations within the lens will produce 
angular broadening, in addition to any broadening resulting from density fluctuations along the 
rest of the line of sight. We shall show that there is a measurable amount of angular broadening 
during the ESE, first by comparing the angular diameter of 1741—038 at different wavelengths 
during the ESE, then by comparing its angular diameter during the ESE to that after the ESE at 
the same wavelength. We shall then relate this additional angular broadening to the fluctuations 
within the lens and show that the lens was probably quite turbulent internally. 

The density fluctuations responsible for interstellar scattering (including angular broadening) 
are typically parameterized by their spatial power spectrum as (Armstrong, Rickett, & 
Spangler 1995) 

Psn^ = Clq~'\ (1) 

The resulting scattering diameter, for a distant source and assuming a = 11/3 (Rickett 1990), is 

Od = 128mas^/(;Hy^SM3/^ 

= I'.'SA^i/^SM^/s, (2) 

where j^gHz is the observing frequency in GHz, Am, is the observing wavelength in meters, and 

SM= j Cl{z)dz. (3) 

In the local ISM Cl ~ lO^^-^ m-^o/s (Armstrong et al. 1995), and SM ~ lO^^-S kpc m^^^/^ ^j. ^ 
typical 1 kpc path length through the local ISM. 

During the ESE we find 13 cm diameters of 0.8 and 1 mas and an 18 cm diameter of 1.6 mas 
(Table 2). Taking the average of the 13 cm diameters, we find the angular diameter scales 
as 9 (X \^'^. We regard this wavelength scaling as consistent with that expected for angular 
broadening, but not decisive evidence in favor of it. We have not attempted to correct for 
any contribution by intrinsic structure, which should have a weaker wavelength dependence 
than angular broadening. Furthermore, some deviation from a strict }?'^ dependence might be 
expected because we are comparing angular diameters determined at different epochs during a 
time-dependent event. Indeed, the angular diameter at 13 cm does change in exactly the manner 
expected if the lens is centrally condensed; near the flux density minimum (1992 June 6 epoch) 
when the path length through a centrally condensed lens would be near the maximum, the angular 
diameter is slightly larger (1 mas) as compared to the angular diameter (0.8 mas) near the end of 



the ESE (1992 July 9) when the path length through the lens would be shorter. However, our 
determination of the A^'® scaling does depend crucially upon the 18 cm diameter, and, in turn, on 
the limited u-v coverage for this observation. If we fit a circular Gaussian instead of an elliptical 
Gaussian to the 18 cm visibility data, we find an angular diameter of 0.3 mas, implying a A~^'^ 
dependence. Hence, we regard the comparison of angular diameters during the ESE as suggestive, 
but not compelling, evidence, for increased angular broadening during the ESE. 

We find a more compelling demonstration of an increase in the angular diameter from 
comparing the angular diameter of 1741—038 during and after the ESE. During the ESE the 13 cm 
diameter of 1741—038 was 0.9 it 0.1 mas. In a series of subsequent observations, during 1994-1997 
(Fey et al. 1996a; Fey et al. 1996b; Fey & Chariot 1997; Fey, unpublished data), acquired while 
the source was not undergoing an ESE, the angular diameter of the compact component has been 
measured to be 0.5 to 0.75 mas, with a mean of 0.63 it 0.04 mas. The quoted uncertainty in 
the mean diameters are statistical. Systematic effects, namely the length of the longest baseline 
in the VLBI array used, probably contribute to an overestimation of the diameter. The same 
systematic effects appear to contribute to our fits also being overestimates of the actual diameter. 
Time-dependent changes may also affect the diameters during the ESE. 

We find the excess angular broadening due to the ESE lens by subtracting in quadrature the 
diameters during and after the ESE. Using the diameters determined above, namely ^in = 0.9 mas 
and ^out = 0.63 mas, we find 5dd = y ^fn ~ ^out ^ 0-''' mas. We treat this value as an upper limit 
based on our assessment of the systematic uncertainties in the measured diameters. However, our 
longest baseline is often significantly longer than that used in determining the diameters after the 
ESE, so our measurement of an increase in the diameter of the source is robust. 

The excess angular broadening of the lens arises from an additional C^ in the line of sight 
to 1741—038 during the ESE. We solve for the level of scattering within the lens, SMiens; in the 
following manner. The scattering within the lens is given by 

SMiens = SMin - SMout, (4) 

where SMin and SMout are the scattering measures seen inside and outside the ESE, respectively. 

We have estimated SMout from RRFID observations (Fey et al. 1996a; Fey et al. 1996b; Fey 
& Chariot 1997; Fey, unpublished data). The RRFID observations acquire 3.6 and 13 cm data 
simultaneously. We have fit the measured 3.6 and 13 cm diameters to 

^ W = (^d,i^m + ^/,l^m> (5) 

assuming that the scattering and intrinsic diameters add in quadrature and that the intrinsic 
diameter scales as A^, as is appropriate for a synchrotron self-absorbed component (Kellerman 
& Owen 1988). Here 9d,i and Oj^i are the scattering and intrinsic diameters, respectively, at the 
fiducial wavelength of 1 m. The spectral index of 1741—038 is a ~ 0.17 {S oc i/") around the time 
of the ESE, but excluding the ESE itself, and a ~ 0.37 near the time of the RRFID observations. 
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in both cases consistent with the assumption of self-absorption. We then use equation (2) to solve 
for SMout- We find SMout = lO"-'^ kpc m-20/3_ 

This value of SMout is consistent with the available constraints on the scattering diameter 
of 1741—038 from low frequencies. Using interplanetary scintillation observations at 92 cm, 
Vijayanarasimha et al. (1985) determined that 1741—038 must have a component with a diameter 
of 100 mas. Assuming that interstellar scattering dominates at this wavelength, we find a 
scattering measure SMout ^ 10~^ kpc m"^'''^. We treat this value as an upper limit because IPS 
observations do not give detailed information on the source structure. 

We estimate SMjn from the difference in the diameter of 1741—038 inside and outside of the 
ESE, 

(se,? = el - Ct = {i'-'8fxT (sMf - sm^/^) . (6) 

We have already constrained 66d < 0.7 mas. We find SMin ~ lO^^'"^ kpc m^^'''^, a value of SM 
that is somewhat larger than the typical SM through the local ISM. 

Thus, the ESE lens toward 1741—038 produced an additional angular broadening of 
approximately 0.7 mas, and the lens itself had a scattering measure of SMiens = 10^^'^ kpc m^^''' ^. 

In order to produce a significant fiux density suppression during the ESE, from stochastic 
broadening alone, the angular diameter of the lens must be comparable to that of the source (Fielder 
et al. 1994a). We assume that the lens diameter was a ~ 1 mas (§3.3) or a ~ 0.1 AU(-D/0.1kpc), 
where D is the distance to the lens. Because the extent of the lens along the line of sight 
may be (considerably) different than its transverse size, we take SMiens = C^ icns'?^' where the 
factor Tj is the ratio of the lens' extent along the line of sight to its transverse extent. We find 
^niens ~ 10^7/"^ m~^'^' ^(D/0.1 kpc)~^. For comparison, Hjellming & Narayan (1986) estimated 
that C^ ^ 10^^'^ ni^^*^' ^ for this line of sight using the refractive scintillation of this source outside 
of the ESE. 

One of the key features of ESE-like events observed toward pulsars is that the pulse width 
does not increase (Cognard et al. 1993; Lestrade et al. 1998; Maitia et al. 1998). Using our 
estimate of SMjens; we can predict how much pulse broadening, another diffractive effect, would 
be produced from a lens comparable to that which passed in front of 1741—038. The amount of 
pulse broadening is (Taylor & Cordes 1993) 

Td < 1.1 ms Z)kpe(SM)6/5z.~22/^ (7) 

We can calculate only an upper limit because the amount of pulse broadening depends upon the 
location of the lens along the line of sight (Cordes & Rickett 1998); the maximum occurs when 
the lens is midway between the observer and pulsar. We predict that ESE lenses typical of the 
one that passed in front of 1741—038 will increase the pulse broadening of a background pulsar 
by only LlD^pc /Us at 1 GHz. This small amount is consistent with the lack of broadening seen 
toward the millisecond pulsars PSR B1937+21 and PSR J1643-1224. 
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Our analysis has assumed that scattering within this lens can be described in terms of a 
power-law spectrum of density fluctuations (eqn. [1]). While this may be true within the lens 
(§3.1.1), the lens itself cannot be formed by the same processes that give rise to the spectrum of 
density fluctuations in the local ISM (Armstrong et al. 1995). First, a medium pervaded by density 
fluctuations on AU scales should produce a correlation between flux density and angular diameter 
(Blandford & Narayan 1985, their Figs. 1 and 4). In contrast, we observe an anti- correlation, with 
the angular diameter increasing as the flux density decreases. Second, if the values we have used 
for a and D are not severe underestimates, the large value of C^ j^^g we infer is significantly larger 
than the value in the local ISM. Even allowing for a structure extremely elongated along the line 
of sight, rj > 100, C„jens remains orders of magnitude above the value in the local ISM. We regard 
the large value of C^ ^^^^ as an indication that the genesis of an ESE lens requires an energetic 
process. 



3.1.1. The Electron Density Power Spectrum within an ESE Lens 
The measured visibility on an interferometer baseline h is 

F(6) = e-^^(^)/2^ (8) 

for a point source seen through a region of density fluctuations with a spatial power spectrum 
given by equation (1). The phase structure function D^{b) is a measure of the fluctuations induced 
in the wavefront's phase as it propagates through the scattering medium and is given by (e.g., 
Cordes & Lazio 1991) 

D4b) c d^y. (9) 

Here 9d is the diffractive scattering angle, and f3 = a — 2. 

There are a number of lines of sight that suggest a ~ 11/3, the Kolmogorov value 
(Rickett 1990). There are also some lines of sight that suggest a significantly larger value, a > 4. 
The diameters tabulated in Table 2 were found by fitting a gaussian to the visibilities, i.e., 
assuming (3 = 2. We now relax that requirement. 

For the three epochs for which we were able to produce images, we fit the visibility data with 
a model of the form of equation (8). For all three epochs we were unable to place any meaningful 
constraints on (3. Allowed values of (3, at all three epochs, were 1-2.2. Though the range is nearly 
centered on the Kolmogorov value of 1.67, we cannot exclude (3 > 2. 



3.2. Refractive Properties of the 1741-038 ESE Lens 

Refractive effects expected from ESE lenses include substructure within the source, angular 
position wander, and multiple imaging, though the extent to which any of these occur depends 
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upon the strength of refraction within the lens. There is no indication of refractively-induced 
substructure in the source such as might be produced if the lens had substructure within it. 
During the event, the source consisted of a compact component, with a weak secondary component 
detectable to the south in at least one epoch. This structure is essentially identical to that seen 
after the event (Shen et al. 1997; Fey et al. 1996a). Below, we quantify the possibility that the 
shape of the lens has altered the shape of the source (§3.4). 

One of the key predictions of a refractive model for ESEs is that an ESE should produce 
angular position wander of the background source. Clegg et al. (1998) predicted that during the 
1741—038 ESE, the angular position of the source wandered by 0.4 mas at 13 cm and by 0.8 mas 
at 18 cm. The observations reported here were not phase-referenced, and the self-calibration we 
performed erased all absolute position information. Furthermore, the proximity of the two intrinsic 
components of 1741—038, the sparse visibility data, and limited dynamic range (particularly 
of the 1992 July 9 observations) give us little confidence of detecting relative position shifts 
between the two components, such as would occur if the lens covered only one component at a 
time. Consequently, even if the lens covered only the brighter component, it would be difficult to 
determine the relative separation, with any degree of confidence, between the two components in 
either the image or visibility domains. We are thus unable to test the prediction of ESE- induced 
angular position wander. 

A second prediction is that an ESE can produce multiple imaging of the background source. 
In the case of this ESE, any secondary image(s) must have been either extremely faint or 
only slightly offset with respect to the primary image, within the 10-15% uncertainties of the 
GBI-measured flux densities (Fiedler et al. 1987b) and those in our amplitude calibration, which 
we estimate to be at least 10%, our models account for all of the flux density measured by the 
GBI. A bright secondary image(s) could have been present only if the multiple images nearly 
overlapped so that no significant anisotropy was produced in our VLBI images (Table 2). Clegg et 
al. (1998) predicted that this ESE was not strong enough to form caustics and produce multiple 
imaging. 



3.3. The ESE Mechanism 

As noted in §1, two general classes of models have been advanced to explain how ESEs 
occur. In the refractive defocussing model (Romani et al. 1987; Clegg et al. 1998; Walker & 
Wardle 1998; hereinafter the RD model) the decrease in the source's flux density during the 
event occurs because of refractive defocussing of the incident rays on the lens. In the stochastic 
broadening model (Fielder et al. 1994a; hereinafter the SB model) the flux density decrease occurs 
because small-scale inhomogeneities in the lens scatter the incident wave front. While most recent 
work has focussed on the RD model — particularly the work of Clegg et al. (1998) who compared 
quantitatively the observed light curve and that predicted from an RD model — both models can 
reproduce the generic features of an ESE light curve, and the only other observational data with 
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which to compare the models have been the pulsar timing programs described in §1. 

A key prediction of the refractive model is that the source's flux density and angular diameter 
should be highly correlated. In contrast we observe an anti- correlation, more consistent with that 
expected from the SB model. Fey et al. (1996b) discussed qualitatively how the SB model could 
produce the 1741—038 ESE, but did not compare quantitatively the observed light curve and 
the SB model. Consequently, we shall re-visit the question of the mechanism by which ESEs are 
produced and whether the SB model can account for this ESE. 

We shall use the SB model developed by Fiedler et al. (1994, Appendix A). This model 
describes the flux density of a source during an ESE as F{t; Fq, fi, 9j, 6i, Of,)- Here Fq is the source's 
nominal flux density outside the lens, n is the proper motion of the lens across the lie of sight, 9i 
is the intrinsic (FWHM) angular diameter of the background source, 9i is the apparent angular 
width of the lens, and radiation incident on the lens is scatter broadened by 6b (FWHM). 

We used a grid-search method to search the available parameter space, evaluating the 
goodness of fit by the x^ statistic. Based on the measurements of 1741—038 outside the ESE 
and our determination of the additional angular broadening during the event (§3.1), we began by 
holding 6 1 and 6^ fixed at 9j = 0.5 mas and Oh = 0.7 mas. We thus fit for the parameters Fq, h, 
and 9i. 

We were unable to find reasonable agreement. In particular the best-fit model had a flux 
density minimum that was approximately 85% of the nominal flux density as opposed to the 50% 
minimum that was observed. Furthermore the shape of the modeled ESE light curve is that of 
a flat-bottomed minimum rather than the rounded minimum observed. We then removed the 
constraints on Oj and 9b separately. If we fit for 9i, Fq, fi, and 9£ while holding 9b fixed, the fit 
agreement improves slightly. However, the depth of the modeled flux density minimum continues 
to be insufficient to match the observed depth, and the minimum is flat-bottomed as opposed to 
the observed rounded minimum. 

If we fit for 9b, Fq, //, and 9i while holding 9i fixed, we find quite reasonable agreement 
with both the depth and shape of the minimum being reproduced. Table 3 lists our best-fitting 
parameters. The significant result of our fit is that the broadening diameter required to reproduce 
the observed light curve is much larger than what we infer from our measurements. The SB 
model alone requires a broadening angle 9b = 2 mas, while our measurements suggest that only an 
additional 0.7 mas of broadening occurred during the event. 

Moreover, the value of 9b found is probably a lower limit. The model of Fiedler et al. (1994) 
assumes that 9b is constant across the width of the lens. If the strength of broadening varies 
across the lens, being stronger in the center and weaker around the edges, an even larger value of 
9b would be required to obtain the same decrease in fiux density during the event. 

We therefore conclude that the SB model alone cannot explain both the observed light curve 
and amount of angular broadening. Given the good agreement that Clegg et al. (1998) found using 
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the RD model, we consider it likely that both refractive defocussing and stochastic broadening are 
occurring. 



3.4. Image Anisotropy 

The images shown in Figures 2-4 display varying degrees of anisotropy. We can identify six 
causes of image anisotropy (Cordes et al. 1987; Cordes 1990; see also Spangler & Cordes 1988): 
(1) incomplete u-v coverage; (2) intrinsic source structure; (3) anisotropic diffractive scattering; 
(4) anisotropic refractive focussing; (5) unresolved multiple images; and (6) a spatially-limited 
scattering volume. 

The first cause — incomplete u-v coverage — is almost certainly the explanation for the image 
of 1992 June 20. The array was elongated significantly in the north-south direction because it 
included Hartebeestock. The u-v coverage was considerably more uniform for the other two 
epochs, though, and incomplete u-v coverage is less likely to be the cause of the image anisotropy 
at these epochs. 

While we cannot rule out the possibility of one or more of the remaining causes contributing 
to the source anisotropy, we can set limits on the shape of the lens by assuming all of the image 
anisotropy is due to refractive focussing. The electron density profile across the lens produces a 
refractive gain G. If the density profile in two orthogonal directions differs, the source will have 
an axial ratio (Cordes et al. 1987; Spangler & Cordes 1988) 

a Ga 

The axial ratios we measure during the ESE are b/a > 0.8. These axial ratios are much closer 
to unity than would be expected if the lens had an axial ratio rj ~ 100 (e.g., Romani et al. 1987), 
and its long axis was in or near the plane of the sky. We conclude that the refractive strength of 
the ESE lens was not considerably different in different directions on the plane of the sky. 

Romani et al. (1987) have speculated that ESE lenses are filamentary structures. Filamentary 
structures have the desirable property of reducing the degree to which the ESE lenses are 
overpressured with respect to the nominal interstellar pressure. As we discuss above, filamentary 
structures would also reduce the inferred level of turbulence within the lenses. Filamentary ionized 
structures would also form naturally in a magnetized medium. If ESE lenses are filamentary 
structures, the nearly isotropic image shapes that we observe indicate that the lenses are extended 
along the line of sight. Lestrade et al. (1998) have suggested that ESEs toward pulsars occur only 
when filamentary or sheet-like structures are favorably oriented along the line of sight. Our VLBI 
images suggest that a similar situation is necessary for ESEs toward extragalactic sources. 
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3.5. Environments of ESE Lenses 

Fiedler et al. (1994b) showed that the hnes of sight to a small number of sources that have 
undergone an ESE display structures suggestive of a turbulent origin. For instance, the line of 
sight to 0954+658 passes near the edge of radio Loop III, with filamentary structures seen in lOO/x 
emission, and the line of sight to 2352+495 (DA 611) passes near RAFGL 5797S, an infrared 
source with a cometary morphology. Fielder et al. (1994a) argued that the distribution of ESE 
sources near the edges of these radio loops was not accidental, but was indicative of a connection 
between ESEs and sites of interstellar turbulence such as old supernova remnants. Supernova 
remnants could also provide a high-pressure environment in which the ESE lenses could survive 
(Romani et al. 1987). Walker & Wardle (1998) have since suggested that ESEs arise from dense 
molecular clouds in the Galaxy's halo. 

The line of sight to 1741—038 displays similar evidence of strong gradients and turbulent-like 
structures. Figure 5a shows the H I column density toward the line of sight of 1741—038 and 
1749+096 (4C 09.57), another source observed to have undergone an ESE. The H I column density 
displays a strong gradient with Galactic latitude, changing by a factor of nearly 10^ over 10° in 
latitude. Furthermore, the two spurs of H I emission extending to higher latitudes in Figure 5a are 
part of radio Loop I, a structure Fiedler et al. (1994a) have already suggested is responsible for the 
ESE toward 1749+096. Like 1749+096, 1741-038 is close to, perhaps within, a portion of Loop I. 
Figure 56 shows the lOO^u emission toward 1741—038. While not as dramatic as the structures 
seen along the lines of sight toward 0954+658 and 2352+495, 1741—038 does appear near a local 
minimum in the 100/U emission. The appearance of the emission is suggestive of a process which 
has excavated a cavity in the ISM. Ionized gas, with significant density enhancements, could result 
within the cavity or on its edges. 

We cannot conclude, from these H I and lOO/j, images alone, that ESEs are produced at 
sites of interstellar turbulence, but these images do add to the existing circumstantial evidence 
suggesting that this is the case. 

It might also be possible to place further constraints on the distance to the material 
responsible for ESEs in this direction. The pulsar PSR J1743-0337 (PSR B1740-03) is located 
only 18' away from 1741—038. It has a dispersion measure of 35 pc cm^'^, corresponding to a 
distance estimate of 1.8 kpc (Taylor &: Cordes 1993). If this pulsar showed enhanced scattering or 
refractive events, like ESEs or fringing in a dynamic spectrum, that would be a strong indication 
that at least a portion of the scattering in this direction occurs because of material closer than 
1.8 kpc. 
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4. Conclusions 

We have presented the first VLBI images of a source, 1741—038, obtained at multiple epochs 
as the source underwent an extreme scattering event. 

We have used these images to assess two models for the origin of ESEs — a refractive 
defocussing model and a stochastic broadening model. The source structure is dominated by a 
compact component and is essentially unchanged during the event as compared to the structure 
after the event. The only change is a slight increase in the diameter of the source (by 0.7 mas), 
an increase we attribute to additional angular broadening within the lens. This additional angular 
broadening is consistent with that expected from a stochastic broadening model but is not 
consistent with that expected from a purely refractive defocussing model. Specifically, a refractive 
defocussing model predicts a correlation between the flux density and angular diameter of the 
source. However, attempts to reproduce the ESE light curve of 1741—038 by a purely SB model 
require a larger increase in the source's diameter (2 mas) than is observed (0.7 mas). We cannot 
test other predictions of the RD model. The refractive defocussing model predicts angular position 
wander of the source, but our observations were not sensitive to absolute angular position shifts. 
We also see no evidence of strong multiple imaging, but a quantitative comparison of the ESE 
light curve and an RD model predicts that the refractive strength of the lens was not sufficient to 
produce multiple imaging (Clegg et al. 1998). We conclude that the 1741—038 ESE involved both 
SB and RD processes. 

The angular diameter of 1741—038 increased by about 0.7 mas during the ESE. The amount 
of angular broadening contributed by the lens implies that the interiors of lenses are highly 
turbulent with levels of scattering orders of magnitude higher than that seen in the local ISM. 
The inferred level of the electron density power spectrum, as parameterized by the coefficient C^, 
is C^ lens '^ lO'^r/"^ m~^'^'^(D/0.1kpc)~"'^. A filamentary lens, with r/ > 1, would decrease the 
required C^. 

The observed visibility data are consistent with the interior of the lens having a power-law 
density power spectrum, with a power spectral index similar to that seen in the local ISM, though 
we cannot rule out a "steep" density spectrum (/? > 4). The lens itself could not have arisen from 
the density fluctuations in the local ISM, however. The value of C^ j^^g is well in excess of the 
local value, and evenly distributed density fluctuations would give rise to a flux density-angular 
diameter correlation in contrast to the observed anti-correlation. 

If ESE lenses are filamentary, as has been suggested to reduce their overpressure relative 
to the ambient medium, ESEs must occur only when the filamentary structures are seen nearly 
lengthwise. A filamentary lens seen transverse to its long axis would produce different refractive 
gains along and across the lens, resulting in image anisotropy. Our images display little anisotropy. 

The line of sight toward 1741—038 shows a strong gradient in the neutral hydrogen density, 
and the source lies close to or within radio Loop I, with small-scale, "cavity" -like structure seen 
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at 100/i. Such interstellar structures along the line of sight to 1741—038 are similar to that seen 
toward some of the other sources for which ESEs have been observed. This line of sight is thus 
consistent with ESE lenses originating from energetic turbulent processes in the ISM. 

Future observations of a source undergoing an ESE will be enhanced by the existence 
of dedicated VLBI arrays such as the Very Long Baseline Array. A key prediction of the 
refractive model — one that we have been unable to test — is the existence of angular wandering. 
Future observations should also have a more extensive and uniform u-v plane coverage, making 
the imaging process easier. Modern, frequency-agile receivers also allow for the possibility of 
simultaneous or nearly simultaneous images at multiple frequencies. Observations at multiple 
frequencies should include imaging the source in the H I line to search for the existence of neutral 
structures related to the ionized structures responsible for ESEs. The major impediment to a set 
of such observations is the lack of a existing monitoring program that could find additional ESEs. 

We thank S. Spangler for helpful discussions regarding interstellar refraction, and K. Desai for 
his image fitting software. This research made use of NASA's Astrophysics Data System Abstract 
Service; the SIMBAD database, operated at the CDS, Strasbourg, France; and NASA's Sky View 
facility^ located at NASA Goddard Space Flight Center. A portion of this work was performed 
while TJWL held a National Research Council-NRL Research Associateship. Basic research in 
radio astronomy at the NRL is supported by the Office of Naval Research. 
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Fig. 1. — The extreme scattering event of 1741—038. The dots show the 2.2 GHz (13 cm) flux 
density as measured by the Green Bank Interferometer. The vertical lines mark the epochs at 
which the VLBI observations reported in this paper were obtained. 



Fig. 2. — The epoch 1992 June 8 at the wavelength 13 cm. (a) The image. The off-source noise 
level is 1 mJybeam"^, and contours are 1 mJy beam~^ x —3, 3, 5, 10, 20, .... The beam is 
shown in the lower left. (6) The visibility data as a function of projected baseline. A model of 
two components is shown superposed for reference. One component is a circular gaussian with 
an amplitude of 1.15 Jy and a diameter equivalent to that implied by the model fit from Table 2, 
namely 1 mas. The second component is a delta function with an amplitude of 0.01 Jy located 
9.7 mas to the south of the first component. 



Fig. 3. — The epoch 1992 June 20 at the wavelength 18 cm. (a) The image. The off-source noise 
level is 1.6 mJy bearn"^, and contours are 1.6 mJy beam^^ x —2, 3, 5, 10, 20, .... The beam is 
shown in the lower left, (b) The visibility data as a function of projected baseline. A model of a 
single, circular gaussian component is shown superposed for reference. The gaussian component 
has an amplitude of 0.95 Jy and a diameter equivalent to that implied by the model fit from Table 2, 
namely 0.61 mas. 



Fig. 4. — The epoch 1992 July 9 at the wavelength 13 cm. (a) The image. The off-source noise level 
is 3.3 mJy beani^^, and contours are 3.3 mJy beam^^ x —3, 3, 5, 10, 20, .... The beam is shown 
in the lower left, (b) The visibility data as a function of projected baseline. A model of a single, 
circular gaussian component is shown superposed for reference. The gaussian component has an 
amplitude of 2.16 Jy and a diameter equivalent to that implied by the model fit from Table 2, 
namely 0.93 mas. 



Fig. 5. — The environs of 1741—038. (a) The gray scale shows the column density of H I (Dickey 
& Lockman 1990), ranging between 10^°'^ cm~^ (white) and 10^^'^ cm~^ (black). The positions 
of 1741—038 and 1749+096, another source observed to have undergone an ESE, are marked, (b) 
The gray scale shows the lOO/U emission, as derived from IRAS observations, ranging between 
17 MJy sr^^ (white) and 37 MJy sr^^ (black). The position of 1741—038 is marked. 
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Table 2. Source Models^ 



Epoch 


A 


S 


r 


i' 


a 


b/a 


^ 




(cm) 


(Jy) 


(mas) 


n 


(mas) 


n 




1992 June 8 


13 


1.15 


0.0 





1.0 


0.83 


-9.5 




0.01 


9.7 


179 


0.0 


1 






1992 June 20 


18 


0.95 


0.0 





1.7 


0.1 


69 


1992 July 9 


13 


2.16 


0.0 





0.98 


0.91 




1992 August 6 


13 


2.52 


0.0 





1.5 


1 





'^Source models consist of gaussians of flux density S, major axis a, and axial ratio b/a at position 
angle (p located a distance r from the phase center at a position angle V'- The stronger component 
was always assumed to be at the phase center. 



Table 3. Parameters of the Best-Fit Stochastic Broadening Model 



Parameter 



Value 



01^ 

Ob 

Fo 



0.5 mas 
2 mas 
1.9 Jy 
1 mas 



14 ;uas d 



-1 



X 



16 



^This parameter held fixed. 
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1.0 
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11 
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1.8 
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3.3 


13 
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f^ ' 1992 June 8 13 0.957 S.5 11 X 4.2 Q -2.2° 1.0 FI 

1992 Juno 20 18 ■ • ■ 6 23 X 7.9 (Q 70° 1.6 Jb 



1992 July 9 13 2.032 6.5 11 X 1.8 <a -4.2° 3.3 FE 

1992 August 6"* 13 2.527 10 ■ ■ ■ ■ • • KI 



'*»»ii,^ , Flux densities obtained from the Green Bank Interferometer monitoring program. Tlie GBI monitored sources only at 3.6 and 13 cm; shown are the flux densities f 

^H ^ closest to the epoch of observation. 

This is the maximum duration for any one baseline, some baselines were on source for shorter periods of time. 
'^Only stations for which correlated data were obtained are listed. Legend: PT, VLBA Pie Town; LA, VLBA Los Alamos; KP, VLBA Kitt Peak; NL, VLBA North Lii 
Hartebeestock; Jb, Jodrell Bank; Mc, Bologna; Nt, Noto; Ro, Deep Space Network DSS65; Ti, Deep Space Network DSS43; Wb, Westerbork 
^^ Data were too limited to map the 1992 August 6 epoch. 



